Scalar self-interactions loosen constraints from fifth force searches 
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The mass of a scalar field mediating a fifth force is tightly constrained by exper- 
iments. We show, however, that adding a quartic self-interaction for such a scalar 



makes most tests much less constraining: the non-linear equation of motion masks 
the coupling of the scalar to matter through the chameleon mechanism. We discuss 
consequences for fifth force experiments. In particular, we find that, with quartic 
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^r) . coupling of order unity, a gravitational strength interaction with matter is allowed 

' by current constraints. We show that our chameleon scalar field results in experi- 

o ' 

mental signatures that could be detected through modest improvements of current 
r-j ■ laboratory set-ups. 

a: 



INTRODUCTION 



We are interested in studying the dynamics of a massive scalar <p with a quartic self- 
interaction and weak couplings to visible matter. The action we wish to consider is 

S= f d'x y/g - -m\<? - - W dsmM) , (1) 

J L J a J -la 

where in the sum over a we give the world-line action for any number of visible particles of 
several species i, whose couplings to the scalar arises from the dependence of the masses rrii 
on 0. Useful dimensionless measures of these couplings are 

, , d log rrii , , 

where Mpi = (87rG) _1 / 2 . The action ([TJ is generally covariant, but in this paper we will 
have occasion only to study effects in flat Minkowski space: g^ u = r]^ u . 

The action ((T} is a standard starting point for discussions of a scalar-mediated fifth force. 
With £ = and to linear order it results in a potential energy between particles of species i 
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FIG. 1: Experimental constraints on the strength a and range A(m) from all fifth force experiments 
to date. Note that a = 1 corresponds to gravitational strength. This plot neglects any self-coupling 
for the scalar field mediating the force. Reprinted from 

and j of the form 

f/(r) = -2A^e^. (3) 

Evidently, the product 2f3if3j is the dimensionless parameter a referred to as the interaction 
strength, while A = m^ 1 is the interaction range. The (3i may be different for different 
species of visible particles, corresponding to isotope-dependence of the scalar forces. Let 
us for now ignore this possibility and take a universal value f3 for all the /3j. Experimental 
constraints on a and A come from a variety of sources. The current state of affairs is shown 
in Figs. H and El [1] . In particular, we see that for a ~ 0(1), corresponding to gravitational 
strength, there is no evidence for a fifth force over the range 100 //m ^ A ^ 10 4 AU. 

If we set £ = 1 in Eq. (0), then the experimental constraints on a and A are drasti- 
cally altered through the chameleon mechanism In particular, we find that the bounds 
from geophysical, LAGEOS satellite, lunar laser ranging (LLR), and planetary data become 
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FIG. 2: Small-scale blow-up of previous figure. Reprinted from P|. 

essentially irrelevant. For instance, the new LAGEOS bound is a < 10 36 . The revised con- 
straints from laboratory experiments are shown in Fig. Our main result is that a ~ 1 is 
allowed on all scales with a quartic self-interaction term. 

For any given m^, it would seem that to allow £ to be of order unity in Eq. ([IJ) is a 
more generic circumstance than constraining it to be negligibly small. Paradoxically, it is 
less generic from the point of view of naturalness in quantum field theory: more precisely, 
Eq. (0) with £ ~ 0(1) and l/m^ ^ 100 /zm is fine-tuned in the presence of otherwise natural 
couplings to the fields of the Standard Model. We will give a more detailed analysis of this 
in section 11111 A failure of naturalness does amount to a serious obstacle to accepting the 
starting point action (0) as theoretically well-motivated. However, we take the view that the 
predictions following from it are sufficiently distinctive and testable that it is worth laying 
them out as an alternative to results based on making <p a free field. Also, in section IIII1 
we explain how the obstacles to naturalness can be to a great extent overcome by positing 
a coupling of (f) to a component of an isospin vector. 

In section |H] we describe the chameleon mechanism as it applies to the model (JTJ and 
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FIG. 3: Revised experimental constraints on the strength a and range A(m) including the effects 
of a quartic self- interaction term with £ = 1. The main difference is for the Eot-Wash [J] and 
Irvine |5|| experiments since their test masses have a thin shell. The test masses used in the 

q n 

Stanford 6] and Colorado [2j experiments, on the other hand, are too small to have a thin shell. 
Thus the corresponding curves are the same as in Fig. [2J 

point out that with £ ~ 0(1), (3 ~ Oil), and = 0, the effective range of the scalar 
force depends strongly on the density of the ambient medium, being roughly a millimeter 
in Earth's atmosphere and roughly 10 km at a generic point in the solar system. This is 
one of the key points which relaxes experimental constraints on a and A when a quartic 
coupling is turned on. Another key point is the failure of the superposition principle and 
the shell theorem that inevitably accompanies the non-linear equation of motion for <p. A 
useful intuition is that for most objects, only a thin shell around the edge of the object 
exerts a significant "pull" on 0, further suppressing the strength of the ^-mediated force. 
In section IIVI we will describe this in more quantitative detail. Then in section IS we will 
describe the modified constraints on a and A with £ = 1. We summarize our results and 
discuss prospects for future experiments in section IV1I 
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II. THE CHAMELEON MECHANISM 

Briefly, the chameleon mechanism is the generation of an effective mass for a scalar 
through the interplay of scalar self-interactions and interactions with ambient matter. Let 
us again assume a universal value j3 for the dimensionless couplings /3, , and let us treat visible 
sector matter as a perfect non-relativistic fluid with negligible pressure. Over-simplifying in 
a way that will not affect the final results, we replace visible sector matter by a single massive 
species with number density n and particle mass m{<p). From Eq. (0) and a hydrodynamic 
approximation one can obtain the following equation of motion for 

™<-, , 9 , £ ,s dm dV e ff 

v eff = + 4i^ 4 + ram (^) 

~ \ m \$ + |j0 4 - j^P + (^-independent) , (5) 

where in the last, approximate equality we have used Eq. (j2J), expanded m(<j)) to linear order 
in (f), and identified nm(<p) as the energy density p. We assume that m(0) is a smooth enough 
function for the linear approximation to be valid for the following discussion. 

Evidently, the equilibrium value of <fi in the presence of non- vanishing p is not 0, but 
rather the minimum <p m i n of V e ff(4>). Small fluctuations around this equilibrium obey a 
Klein-Gordon equation with a mass given by 

2 _ d 2 V eff ((f) min ) _ _ 2 g 2 ^ d 2 m((f) min ) 

where the last term would be neglected in the linear approximation for m(<ft). Clearly, the 
effective mass is larger for larger p: the scalar force is screened, because of the quartic 
self-interaction, in the presence of matter. 

This chameleon mechanism is most spectacular when m^ = 0. Then, using the linear 
approximation for m(<j)), 

Q(3p n 1/3 



M Pl (7) 
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In more user-friendly units: 

1/3 



■'linn 



a" 1 ] « 10 (0 (pfg/cm 3 ]) 1 / 3 (9) 
m^H « 0.1 (^ 2 )- 1/6 (p[g/cm 3 ])- 1/3 . (10) 

Corrections to these formulas would be suppressed by factors of pjMp X or mi „ <C M Ph both 
of which are indeed small quantities for any reasonable matter density. 

For example, the atmosphere has mean density p « 10~ 3 g/cm 3 . Substituting in Eq. (fTUj) 
and assuming /3,£ ~ C(1)j we find 



/;/ 



-l 

atm 



1 mm . (11) 



Hence, even though the field is exactly massless in vacuum, we see that the force it mediates 
has a range of one millimeter in the atmosphere. 

In the solar system, the relevant matter background is the nearly homogeneous baryonic 
gas and dark matter, with density pc ~ 10~ 24 g/cm 3 . The corresponding interaction range 
is then 

m^ 1 w 10 km. (12) 

In jsj], the chameleon mechanism relied upon the interplay of a runaway potential for <fi 
and exponential couplings to matter, and one scale had to be set by hand in order to achieve 
something similar to Eq. (|llj). Thus, a novelty of the present setup is that millimeter range 
screening comes to us "for free" given a quartic self-interaction with a coefficient of order 
unity. Another distinction in comparing with jjj is that the ^-mediated interactions are 
short range both in the atmosphere and in the solar system. 



III. NATURALNESS OF THE c/> 4 CHAMELEON CONSTRUCTION 



Consider the following alternative to our starting point action (jTj): 

"l,~.,o 1 



S 



d x 



2 ±2 



— 4 + ipi$ip — m((f))ipip 



(13) 



where ip is a Dirac fermion. Once again, we have over-simplified considerably by replacing 
visible sector matter with a single free fermion, but the action ()13)) will nevertheless serve 
to illustrate points that carry over to a more general setting. At the classical level, for 
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non-relativistic processes at low enough density that the Pauli Exclusion Principle is not 
important, Eq. (JT3|l is equivalent to Eq. (JTJ). It is assumed that m(0) = m{Q)h(f3(f)/Mpi) for 
some function h(x) = 1 — x + 0(x 2 ), whose further Taylor coefficients are at most of order 
unity. 

If £ = = m^, then integrating out ip via graphs of the form shown in Fig. leads to 
a potential for of the general form V e ff(4>) = A 4 g(f3(j)/Mpi) for some function g(x) whose 
values and derivatives are of order unity (apart from logarithmic singularities at points 
where m(0) = 0). A local minimum of Ks(0) would typically give a mass for of the form 
~ A 4 /3 2 /M Pl . If A ~ 5 TeV — a reasonable scale for present purposes since we have no 
direct knowledge of physics above the TeV scale — then m e //[meV] ~ 10/3. (The coefficient is 
slightly lower if one accounts for phase space factors |2j.) This is a standard line of argument, 
whose conclusion is that a range on the order of a millimeter for a scalar-mediated force is 
natural given a cutoff near the electroweak scale and Planck-suppressed couplings. 

Modern-day experiments probe a much broader range of lengths, despite the argument 
just summarized that naturalness favors the sub-millimeter regime. Let us then regard 
naturalness as a guide to favored regimes but not as a firm rule. In this spirit, it is instructive 
to see how naturalness becomes more difficult to arrange in a theory with £ of order unity. 
At the end of this section, we will suggest a way to circumvent the main difficulty. 

As a warmup, let's first consider the case where £ ~ 1/10 but the fermions are decoupled 
(i.e., dm/d(f) = 0). Then having small or zero is no less natural than when £ is 10]. 
This seems to conflict with the prediction that the self-energy graph in Fig. |3Jd contributes 
£A 2 to m|. Let us therefore analyze the issue carefully in a renormalization group (RG) 
language. 

Assume then that £ = £(A) = 1/10 and = in the action that best describes the 
dynamics at a scale A (this is an effective action in the Wilsonian sense, and the specific 
value of £(A) is chosen so as to have some control in perturbation theory). Working to 
one-loop order, the effective potential (part of the 1PI effective action) is (see for instance 
) 

K//(0) = ~rr<P + TT7r-r 2 <\> Hog -J - -r , (14) 



4! r (16vr) 2 ^ \ [i 2 6 
where the renormalization conditions are ^pr(O) = and ^r(/i) = £(/-0, and the fact that 
V e ff is independent of the arbitrary scale /i is encoded in the RG equation = j§^i£, 2 - 
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(a) (b) 

FIG. 4: Loop diagrams for the action (|13|) which contribute to the scalar potential. 

A difficulty with (fT4"j) is that the 4 log <fi 2 term causes symmetry breaking: the minimum 
is at some non-zero 0, and <fr at this minimum has a mass. This is easily cured by adding a 
fermion, call it \i with mass m x (<f)) = r](p, where r\ is some additional coupling of order unity. 
This fermion is not part of the Standard Model: it is an extra field whose sole purpose is 
to reverse the sign on the 4 log0 2 term to insure a stable minimum at <fi = 0. The main 
point, which the fermions do not alter, is that is absent at lower scales if it is absent at 
a high scale. If ^ 0, then it evolves according to an RG equation: 

^ = (li^ + 4^) m ^ (15) 

This multiplicative renormalization of may seem in conflict with the idea that the self- 
energy graph in Fig.HjD tends to "drag" up to the cutoff A. The resolution is that adding 
£A 2 to is (part of) how one gets from a classical action to a quantum effective action, 
whereas the renormalization equation (|15jl is how one gets from a Wilsonian action for the 
quantum theory at one scale to the effective action at a lower scale. 

From the Wilsonian perspective, then, before adding fermions, non-zero £ with = 
is indeed more generic than £ = = m^. Super symmetry and / or string theory however 
provide numerous examples of supersymmetric field theories with flat directions or string 
moduli, for which £ = = may indeed be nearly right at some high energy scale. How 
the concept of genericity may be altered in string theory is only beginning to be understood. 

There are examples in supersymmetric field theory where the renormalizable scalar po- 
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tential is purely quartic: for example, consider a U(l) gauge theory coupled to three chiral 
superfields, one of charge 1 and the others of charge —1/2. Gauge invariance restricts the 
superpotential to purely cubic terms. The potential receives quartic contributions from 
D-terms and F-terms, and cases with a unique global minimum at the origin are readily 
constructed. It is more difficult to couple such a theory to the MSSM in such a way that the 
masses of quarks m q (<j)) give (3 = — Mpid log m/dip on the order of unity, where denotes 
the scalars charged under the extra U(l). 

Let us put supersymmetry aside and simply assume that the scalar potential is |^0 4 at 
some high scale A. The real problem arises when we include loops of Standard Model fields, 
because the latter tend to change the potential in such a way as to shift the vacuum away 
from = 0. To see this, it is enough to go back to the action and compute the one- loop 
effective potential: before any renormalization, 

V eff {4>) = W0) + 3^ B- 1 )^™^) 2 - £(-l)*W) 4 log > ( 16 ) 

where the sum is over the degrees of freedom (counting a real boson as 1 and a complex 
Dirac fermion as 4), and (— 1) F is 1 for a boson and —1 for a fermion. The quantity 
appearing in (fTp^l is by definition d 2 V e ff(<j))/d<j) 2 . Evidently, V e ff(4>) has odd terms in even 
if Vbarei.^) doesn't. So = isn't even an extremum of V e ff((p), and at the true extremum 
will have a mass as well as cubic and and higher interactions. We do not see how to make 
an interesting theory of chameleon fields in this way. 

The problem originates in the loss of the Z 2 symmetry — > — that arises when we 
couple to a fermion with mass m(0) = mo(l + /30/Mp;). A chiral phase rotation ip — > e l7r75, 
can switch the sign of mo but not of /3. (This explains why the fermion x mentioned earlier, 
with m x ((j)) = r](p, doesn't break the Z 2 : one just extends its action to include x ~^ elnl5 X-) 
Let us try to get closer to a field-theoretically natural chameleon by considering coupled 
to up and down quarks through a term that is part of an isospin triplet: 

S = J d^x^fg Q(<90) 2 - 7^j0 2 - |j0 4 + x{i0 ~ V<t>)x + ~ m )l ~ fyqnq^J , (17) 

where q denotes the isospin quark doublet (^), and r 3 = ( _ x ) is one of the Pauli matrices. 
The action ()17)) is still conspicuously imperfect: the quarks are free except for the interaction 
with 0, and we have given them an equal mass m. But the <fi — > —0 symmetry can now 
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be extended to include an isospin rotation sending u — *> d and d — > —u, and the effective 
potential, for ^/rj 2 not too large, is even with a minimum at <fi = 0. 

In a more realistic treatment including QCD, there is chiral symmetry breaking through 
an isospin-singlet VEV (qq), but (qrq) = to avoid breaking isospin. The absence of an 
isospin triplet VEV is fortunate because otherwise there would be a tadpole for <j) that would 
lead us back to V e ff(4>) that is not even. 

If isospin were a perfect symmetry of QCD, then the coupling —/3<f)qT 3 q would result 
in a naturally even scalar potential with small preserved by RG flow, and opposite 
dimensionless couplings $ (defined as in (J2J)) for the proton and neutron, proportional to (3. 
Isospin is slightly broken by explicit quark masses and electromagnetic effects, and we defer 
to the future an attempt to determine how much this alters the final picture. 

In summary, the simplest model that we study, where = and the j3i for visible 
particles are all equal to a common value, is not field-theoretically natural because of terms 
in V e ff(4>) that violate <p — > — symmetry. One could contrive to eliminate such terms by 
cancelling contributions to V e ff(<f)) from the first generation of quarks against contributions 
from the other generations: this is obviously a fine-tuned approach since the interactions of 
the heavy quarks are quite different from u and d. An alternative model, with equal and 
opposite f3i for protons and neutrons, is more promising from the perspective of naturalness. 



IV. REVIEW OF THE THIN-SHELL MECHANISM 

In the previous sections we have encountered one of the characteristic features of 
chameleons, namely that their mass depends on the surrounding density. Another important 
property of such scalars is that their effective coupling to matter is also density- dependent. 
To be precise, the 0-field outside a body involves an effective coupling, (3 e ff, which, for suffi- 
ciently dense objects, is much smaller than the "bare" coupling (3. Thus the 0-force between 
two dense bodies is suppressed by this effect, first referred to as a "thin-shell" mechanism 

ings between 



0. 



in P| for reasons that will soon be obvious. Density-dependent effective coup 
a scalar field and matter were first observed by Damour and Esposito-Farese 

Consider a spherical body of mass M c , radius R c and homogeneous density p c . The body 
is assumed immersed in a homogeneous medium of density poo. Throughout we denote by 
and the field value that minimizes V e ff for p = p c and p^, respectively. Similarly, the 
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mass of small fluctuations about these minima are m c and m^. Equation Q thus reduces 
to 

^ + r^ = m ^ + 3! ~M^ P(r) ' (18) 
where p = p c for r < R c and = for r > R c . This differential equation is subject to the 
following boundary conditions: i) the solution must be regular at the origin; ii) far away 
from the body, the field must tend to its expectation value in the ambient medium. The 
latter ensures that the force on a test particle becomes vanishingly small as the test particle 
moves infinitely far. Thus, we impose 

^ = at r = ; 
dr 

— > 0oo as r — > oo . (19) 

Although the above equation is non-linear, we can use the linear approximation which, 
as we will see, is valid for all objects of interest. For r < R c , substituting = C + 5(f) in 
Eq. (HU gives 

d 2 5(t) 2 dScj) 



i • +-~r = ™M- (20) 
dr z r dr 

The solution satisfying the first boundary condition above is 

„ . Asinhm c r , . . 

0(r < R c ) = — ^- + <p c , (21) 

where A is a constant to be determined later. 

Similarly, far away from the body, we can once again linearize Eq. (|18|) . this time about 
0oo. The solution satisfying (f> — > (j)^ as r — > oo is 

0(r > R c ) = + , (22) 

where B is another constant. 

Strictly speaking, this solution is only valid in the asymptotic region far from the object. 
Nevertheless, to obtain an approximate analytical solution, we will assume that it holds 
even near the surface, r ~ R c . We will later determine the degree of accuracy of this 
approximation by comparing with numerical calculations. Imposing continuity of <fi and 
dcj)/dr at r = R c then allows one to solve for A and B. This gives 

0(r < R c ) = I \+ m ~ R c X We - <U smhKr) + ^ . (23) 



m c R c coth(m c i? c ) + m^Rf, J sinh(m c i? l 
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0(r > R c ) = { ^RcCoMmM - 1 1 Wc - M ^r-y + ^ (24) 



m c i? c coth(m c i? c ) + m^Rc 
This, however, assumes linearity for r < i? c . It is easily seen from Eq. (123(1 that this is a 
valid assumption provided that 

m CJ R c > 1 . (25) 

If we further assume that the density contrast between the object and the ambient matter 
is high, so that (p c ^> (p^ and m c ^> moo, the exterior solution in Eq. reduces to 



ft tf M p-^oo(r-iio) 
0(r > i? c ) « + 0OO , (26) 

with 

where in the last step we have used Eqs. (jZJ and ( |5|l . Equation is precisely the exterior 
solution for a normal scalar of mass moo, with /3 e // replaced by (3. Thus, due to the quartic 
self-coupling, the exterior solution for objects satisfying m c R c 3> 1 is exactly that a point 
particle albeit with a much weaker effective coupling constant, (3 e ff <C (3. The physical 
interpretation is clear: non-linear interactions responsible for this suppressed coupling are 
only effective provided the object is sufficiently large compared to the Compton wavelength 
of the chameleon, so that the latter can essentially "feel" the object. 

n 

This suppression mechanism was called "thin-shell" in |3|. This terminology is motivated 
by the realization from Eq. (|21j) that (ft stays essentially constant throughout the bulk of the 
object except within a thin shell of thickness ~ m" 1 . Thus only this thin shell contributes 
to the force on a test particle, resulting in a suppressed effective coupling /3 e ff- 

To corroborate these statements, consider the opposite limit, m c R c < 1. In this case, 
the linear approximation does not hold, and therefore <ft *C <f) c for r < R c - Equation (fTHJ) in 
the range < r < R c thus becomes 

d 2 (f) 2# p 

with solution 

do T 2 

<f>(r < R c ) = -j£r- + const. (29) 
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FIG. 5: Profile for the chameleon for an object with m c R c ~ 65. The solid curve is the result of 
numerical integration; the dotted line is the analytical approximation derived in the text. 



The exterior solution (r > R c ) is once again of the form 0(r) ~ e 
dcj)/dr at the surface of the object gives 



(r > Rr) « 



+ 



"jr. Matching and 



(30) 



AnM Pl r 

As anticipated, the exterior field is that of a point particle with unsuppressed coupling in 
this case. Moreover, there is no thin shell since, as seen from Eq. ()29|). the field never remains 
nearly constant within the object. 

The above analytic expressions have been checked against numerical calculations. Con- 
sider P = 1, £ = 1, p c = 1 g/cm 3 and = 10~ 6 g/cm 3 , corresponding roughly to a ball 
of Beryllium in thin air. Correspondingly, from Eqs. (fTUj) we have <p c ~ m c w 10 mm -1 
and 0oo ~ moo ~ 0.1 mm -1 . Figure El shows the result of numerically integrating Eq. ()18|) 
for an object of radius R c ~ 6.5 mm. Since m c R c ~ 65, we expect a suppressed coupling 
in this case. This is confirmed by the numerics. The solid curve in Fig. El is the numerical 
solution. We see that indeed the field remains at <fi ~ C for < r ^ R c , justifying the linear 
approximation in this regime and supporting the "thin-shell" interpretation. The dotted line 
is a plot of Eqs. (J23*|) and ()24|). For r > R c , the exact solution is initially steeper than the 
analytic approximation; very quickly, however, the roles are reversed, and it is the dotted 
line that it is steeper. The two curves eventually both decay exponentially. 

The small discrepancy for r ^ R c is a consequence of our trusting Eq. near the 
surface of the object. To assess how this affects our predictions for the suppressed coupling, 
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FIG. 6: Ratio of the numerical to theoretical chameleon-mediated force on a test particle in the 
exterior of the object. 



let us define 



n{o,ctual) n 



^ ^(actual) 
£ff \7(f)(theo) 



(31) 



where the superscripts (actual) and (theo) refer to the numerical and analytical solutions, 
respectively, and where /3 e // is given in Eq. (|27|). Thus (3 { e a f c } ual) /{3 eff is the ratio of the 
numerical to theoretical mediated force on a test particle. Fig. |U]is a plot of this ratio. We 
see that the discrepancy ranges from « 7 at r = R c to « 0.1 for r 3> -R c - In other words, the 
suppressed coupling is found numerically to lie in the range 2 • 10~ 4 ^ ^{actual) < 1 5 . 1CT 2 , 
while the analytical estimate in Eq. (|2*7)l gives /3 e // ~ 2 • 10~ 3 . The discrepancy disappears, 
of course, as m c R c — > 1. In particular, the force between two macroscopic test masses is 
always less than or equal to the force between two point particles of the same mass. 



V. MODIFIED FIFTH FORCE CONSTRAINTS 



We have seen that a quartic interaction suppresses the 0-mediated force in two ways. 
Firstly, the ambient matter density gives a non-zero expectation value to about which the 
mass of fluctuations can be much larger than in vacuum. Secondly, non-linear interactions 
of (j> inside a dense macroscopic body leads to a suppressed effective coupling as felt by an 
outside test mass. In this section, we explore the consequences of these two effects for fifth 
force searches. 

Consider two bodies of radius Ri and mass Mi, i = 1,2, both satisfying Eq. (J23J). For 
simplicity, we assume that they have equal bare coupling (3 and density p = 1 g/cm 3 . The 
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latter is the ballpark value for all bodies used for fifth force constraints, let it be test masses 
in the laboratory, satellites or planets. From the second of Eqs. (jlOj). this corresponds to 
m~ l « 0.1 (£/3 2 ) _1//6 mm. The potential energy associated with the fifth force is then 

U(r) = , ,32) 

where a e ff is given by Eq. (|27)l and the relation a = 2f3 2 \ 

-2 ( a \ 1 

a ?ff ~ 10 72" "B27 T52T T ' ( 33 ) 

\£ / -Rf[mmJi?2[ mm J 

and where X e ff is the effective Compton wavelength of </> in the surrounding medium. From 

Eqs. (JBJ) and f|TTT|> - we find for the latter 

\ eff = \ . (34) 

1 + 100(C/3 2 ) x /3p^ [g/cm 3 ]A 2 [mm] 

A. Laboratory Searches 

Laboratory fifth force experiments are performed in vacuum, m 0, corresponding to 
an unsuppressed interaction range: 

Kff = A . (35) 

The tightest constraints from the laboratory on forces mediated by scalars with insignifi- 
cant self-interaction come from Eot-Wash |4| and Irvine [5J. Typical test masses for these 
experiments have characteristic size Ri ~ 10 mm, resulting in 

1/3 

?2 j • (36) 

To be completely accurate, however, in laboratory experiments the separation between test 
masses is often smaller than their size. We have seen in Sec. II VI that, in this regime, the 
prediction for (3 e ff is smaller by a factor of 10 or so compared to the numerical solution. 
Including this fudge factor, we obtain 

1/3 

e ) • (37, 

It is the effective coupling a e ff that is truly constrained by experiment. Equation (|3Tj) thus 
implies that the actual coupling strength a can be much larger. With this revised perspective 



_ 4 / a 
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on the data, Fig. ^ for instance, says that laboratory experiments constrain a e ff ^ 5 • 
(from Irvine e 5]), implying the much weaker constraint on the bare coupling a ^ 10 2 for 
£ ~ Oil). In particular, we should stress that a ~ 0(1), the expected outcome from string 
compactifications, is allowed. 



B. Planetary, Lunar and Geophysical Constraints 

The existence of a fifth force is also constrained by solar system data, in particular from 
the LAGEOS satellite, Lunar Laser Ranging (LLR), and anomalous perihelion precession 
of planets. See Fig. ^ For all of these, the effective interaction range is determined by the 
ambient baryonic gas and dark matter, with approximate density pa ~ 10~ 24 g/cm 3 . Thus 
Eq. (jmj) implies 

A e// <10(63 2 )- 1/6 km. (38) 

Therefore, even for £ ~ 10~ 6 (and reasonable (3), this gives A e // ^ 100 km. From Fig. d 
on the other hand, we see that LLR and planetary data only significantly constrain a fifth 
force for ranges roughly greater than 1000 km. Hence, unless £ is exceedingly small, the 
constraints from LLR and planetary orbits become completely irrelevant. 

One is therefore left with the LAGEOS data, a satellite orbiting the Earth at an average 
distance of approximately 10 4 km. The radius of the Earth is of order 10 10 mm, while that 
of the satellite is ~ 10 3 mm, resulting in an effective coupling strength of 

«.// - 10- 28 V3 . (39) 

We infer from Fig. ^ that the LAGEOS constrains a e // ^ 10~ 4 for A e /f ^ 100 km. The 
above thus implies that this bound is consistent with a bare coupling constant as large as 
a ~ 10 72 e 2 ! 

Geophysical constraints result from gravity experiments performed in mines, boreholes, 
lakes, oceans and towers. All of these essentially measure the relative difference in the 
gravitational force on two test masses at different depth. Since this involves the gravitational 
field of the Earth, substituting R & ~ 10 10 mm in Eq. (|3^j) gives 

« e //<10- 22 (|r) 1/3 . (40) 

Once again it is clear from Fig. ^ that these constraints can thus be neglected for reasonable 
parameter values. 
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C. Isospin triplet coupling 

In section IIII1 we concluded that the difficulties in constructing a model that satisfied 
field theoretic notions of naturalness could be ameliorated if (3 were equal and opposite for 
protons and neutrons (and zero for electrons). Such a situation can be termed an isospin 
triplet coupling of the scalar to matter. Let us briefly review how this model compares with 
our main focus, namely a universal value of j3 for all particles. 

Clearly, the effective screening length m~jf now depends not on the ambient density p, 
but rather on the difference p v — p n of the density p p due to protons and p n due to neutrons. 
The key result (flT)J) carries over with the replacement p — > \p p — p n \. In Table H] we list some 
representative values of rn~^ and /3 e // for an object 1 cm in radius, assuming £ = 1/10 and 
(3 = 1. Dry air has an excess of neutrons, due mostly to the fact that argon makes up about 



Material 


p [g/cm 3 ] 


\Pp ~ Pn\/P 


m eff 


f3 e ff/(3 for R c = 1 cm 


lithium 


0.4 


0.14 


0.3 mm 


8 • 10"3 


aluminum 


2.7 


0.04 


0.3 mm 


8 • 10-3 


copper 
gold- 


8.9 
19.3 


0.1 
0.2 


0.15 mm 
0.1 mm 


2 • 10-3 
9 • 10" 4 


dry air 

interstellar 
medium 


1.3 x lO" 3 

1Q -24 


1.3 x 10^3 
1 


1 cm 
10 km 





TABLE I: Effective interaction range and coupling in various materials for the isospin triplet model 

a percent of air and has a 10% excess of neutrons over protons. Humid air picks up extra 
protons from water molecules. This results in \p p — p n \/p = when the dew point is roughly 
17 °C, corresponding to a relative humidity of approximately 80% at 20 °C. 

In light of the sub-millimeter screening lengths in the first several entries of Table HI 
laboratory experiments of the type described in section IV Al result in constraints that are 
not much different from the case where (3 is the same for all particles. Indeed, most test 
masses used in fifth force searches are made of materials such as copper, beryllium, iron etc., 
all of which have a proton overdensity of about 10% as seen from Table HI This translates 
into an increase of about 10 2//3 « 4.6 in the effective coupling /3 e //> as seen from Eqs. (jHJ 
and (j2"Tj) . All the scaling laws, e.g., a e ff ~ (a/£ 2 ) 1//3 as in Eq. (J27J), still hold. In particular, 
a ~ 1 is still allowed, provided that £ <: 10 2 . Meanwhile, the analysis of planetary, lunar 
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and geophysical constraints of section IV Bl is only mildly altered since the Earth consists 
primarily of iron with proton overdensity of about 5%. Thus, a e ff remains ridiculously 
small, as in Eqs. (jSHJ and 



VI. SUMMARY AND DISCUSSION 

Figure El shows the revised constraints on a and A for £ = 1. The key result is that 
a ~ 0(1), corresponding to a fifth force of gravitational strength, is allowed for all scales 
probed by experiments. Our model essentially escapes the Eot-Wash and Irvine bounds 
because the test masses used in these experiments are large enough to suffer from a thin- 
shell suppression. 

Nevertheless, we wish to stress that a modest improvement in current laboratory set-ups 
could lead to an unambiguous detection of the chameleon. Indeed, consider a laboratory 
experiment to detect a fifth force between two test masses of radius R in vacuum. We have 
in mind a set-up similar to the Eot-Wash experiment, for instance. Assuming (3 = 1 and 
£ 1, then Eq. (J33|) gives 

a eff ~ Kr 2 iT 4 [mm] . (41) 

The Eot-Wash experiment used test masses of characteristic size R ~ 1 cm. This gives 
a e ff ~ 10~ 4 (see Eq. (j37j) ). which is just under the radar of current experiments, as detailed 
in Sec. IV Al Repeating the same experiment with test masses of characteristic size R ~ 1 mm 
would instead give 

a e// ~l(T 2 , (42) 



which is well within the sensitivity range of the old 



lot- Wash experiment. 



Experiments such as Stanford Jf| and Colorado 3| have used test masses smaller than 
1 mm, which therefore do not suffer from any thin-shell suppression. Indeed, comparing 
Figs. El and El we see that the corresponding curves are left unchanged by the addition of 
a quartic self-coupling. In particular, Fig. El shows that the Colorado experiment provides 
the tightest bound on a on small scales. Modest improvements in its sensitivity would allow 
Colorado to probe the theoretically interesting range of a ^ 1. 

Thus, by making the test masses a bit smaller or improving the sensitivity, one might 
detect the fifth force described here — provided the couplings (3 are not significantly smaller 
than the natural value, namely unity. This detection would unambiguously be identified as 
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resulting from a chameleon field, since the above fifth force would only appear for sufficiently 
small test masses. This could easily be verified explicitly, by performing the experiment with 
different size test masses and checking for the R~ 4 dependence predicted above. Moreover, 
as shown in Fig. the fifth force is steeper than 1/r 2 at short distances. (If the force 
were 1/r 2 at short distances, then p( actual ) j p would be constant.) This should result in a 
torsion-pendulum signal that is distinguishable from the usual Yukawa fifth force. 

Because it is possible in the near future to start probing the interesting range of param- 
eters for our model, it would certainly be helpful to refine our calculation of the expected 
chameleon force and make it more tailored to a realistic experiment. For instance, test 
masses used in an Eot- Wash-like experiment are not spherically symmetric, as assumed 
here, but instead cylindrical slabs. We leave for future work the calculation of the precise 
force and expected signal for such a set-up. 
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